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The energy dependence of the total hadronic cross sections is caused by gluon bremsstrahlung which 
we treat nonperturbatively. It is located at small transverse distances about 0.3 fm from the valence 
quarks. The cross section of gluon radiation is predicted to exponentiate and rise with energy as s A 
with A = 0.17 ± 0.01. The total cross section also includes a large energy independent Born term 
which corresponds to no gluon radiation. The calculated total cross section and the slope of elastic 
scattering are in good agreement with the data. 
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The dynamics of energy dependence of the hadronic to- 
tal cross sections is a long standing challenge since 1973 
when this effect was first observed at the ISR. In DIS 
the source for the rising total cross section for interac- 
tion of highly virtual photons is well understood in QCD 
as caused by an intensive gluon bremsstrahlung . In- 
deed, radiation of each gluon supply an extra Ins and 
InQ 2 . This is a specific regime of radiation when a qq 
fluctuation of the photon of a tiny size ~ 1/Q radiates 
gluons at much larger transverse separations. 

It is difficult to extend the perturbative results to soft 
hadronic collisions because it is quite a different regime 
where the approximations made in the perturbative case 
break down. Namely, gluon radiation giving rise to 
the energy dependence of the total cross section occurs 
at rather small transverse distances around the valence 
quarks, tq ~ 0.3 fm which are much smaller than the 
mean interquark spacing in light hadrons. This conclu- 
sion follows from the analysis Q of the data for diffractive 
gluon radiation based on the light-cone approach when 
the effective nonperturbative interaction of radiated glu- 
ons is included. 

The smallness of the gluon clouds of the valence quarks 
is confirmed by the study of the gluon formfactor of the 
proton employing QCD sum rules 0. The Q 2 depen- 
dence of the formfactor turns out to be rather weak cor- 
responding to a small radius of the gluon distribution 
which was estimated at the same value rg « 0.3 fm. 

Another evidence for a short gluon-gluon correlation 
length A ~ 0.3 fm arises in the stochastic vacuum model 
of Dosch and Simonov |||| , as it was measured on the 
lattice jjj . In the case of a poorly populated gluon cloud 
(only about one gluon is radiated by a valence quark 
at available energies, see below) this corresponds to the 
correlation radius between the quark and the gluon. 

The same size ~ 0.3 fm emerges from the Shuryak's 
instanton liquid model g as the instanton size which 
controls the mean radius of the sea surrounding a valence 
quark and by many phenomenological analyses. 

In the Gribov's theory of confinement P,|l0fl the same 



distance ~ 0.3 fm should correspond to the critical 
regime related to breaking of chiral symmetry. Namely, 
at smaller distances, a perturbative quark-gluon basis is 
appropriate, while at larger separations quasi-Goldstone 
pions emerge. The corresponding critical value of the 
QCD running constant a c = 0.43 evaluated in Q turns 
out to be very close to our estimate (see below) of a s cor- 
responding to gluon radiation separated by 0.3 fm. The 
value of a s is crucial for our evaluation of the energy 
dependence of the gluon bremsstrahlung. 

It is quite plausible that all these observations are the 
manifestations of the same dynamics, however it is still 
unclear how to make a Lorentz boost in these approaches. 
This is the advantage of the light cone treatment of non- 
perturbative gluon radiation || which seems to be best 
designed for calculating the energy dependence of the to- 
tal cross section. We believe that the nonperturbative 
interaction of gluons introduced in ||] as a light-cone po- 
tential is an effective manifestation of properties of the 
QCD vacuum. Similar scale ~ 0.3 fm found in all these 
approaches supports this conjecture. 

An interesting attempt to implement the nonpertur- 
bative gluon interaction into the Pomeron ladder build- 
ing was made recently by Kharzeev and Levin Jlifl and 
Shuryak jl5| . They found that the radiation of colorless 
pairs of gluons is a part of the leading-log approximation 
since each extra power of the coupling a s cancels due 
to the strong glue-glue interaction. The radiated glue- 
balls are not clustering around the valence quarks, but 
spreading all over the hadron. The estimated A ss 0.05 
p5]| is about twice as small (and even more so if corrected 
for unitarity) as the data need. Although the scale for 
a' p ~ V-Mq seems to be correct, an extra factor A/4 
makes it too small. 

We start calculating the energy dependence of the total 
cross section summing up the contributions of different 
Fock components of the incident hadron, 
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To avoid double-counting, we sum over cross sections a n 
of physical processes corresponding to the radiation of n 
gluons. 

The lowest Fock component of a hadron contains only 
valence quarks. The corresponding Born term in the total 
cross section has the form (for the sake of simplicity we 
assume that the incident hadron is a meson), 
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Here the Fock state wave function ^/ qq (a q , R) depends 
on the transverse q — q separation R and on the fraction 
a q of the light-cone momentum of the pair carried by 
the quark. The cross section a qq (R) of interaction of the 
valence qq dipole with a nucleon cannot be calculated 
pcrturbatively since the separation R is large. According 
to | fL3f this energy independent term has no relation to 
the smallness of the spots (gluon clouds) in the hadron. 

The next contribution to comes from the radiation 
of a single gluon. The radiation is possible only due to 
the difference between the cross sections for the qq and 
qqG Fock components, otherwise no new state can be 
produced ||. The cross section of radiation of a single 
gluon reads S, 
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Here a G is the fraction of the quark momentum carried 
by the gluon, and r is the quark-gluon transverse separa- 
tion. The three terms in the curly brackets correspond to 
the radiation of the gluon by the quark, by the antiquark 
and to their interference respectively. 

The nonperturbative wave function for a quark-gluon 
Fock component is derived in ||. Neglecting the quark 
mass, the wave function reads, 
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where e is the polarization vector of the massless gluon. 
The parameter &o = 0.65 GeV characterizing the non- 
perturbative quark gluon interaction is fixed by the data 
on large mass diffractive dissociation corresponding to 
the triple-Pomeron limit. It leads to quite a short mean 
quark-gluon separation tq = \J (r 2 ) = l/&o ~ 0.3 /m, 
which is small relative to the hadronic size. Therefore, 



only one or the other of the first two terms in ([|) can be 
large, while the interference one can always be neglected. 
In this case, the integration in (||) is easily performed, 
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Here we assume that the approximation <Jq q {r) — Cr 2 is 
valid for r ~ l/^o- N is the number of valence quarks, 
m(s/s ) = ln[(a G ) mal /(a G ) mm ], where (a G ) min = 
2 6§/s, but (a G ) max is ill defined. It should be suffi- 
ciently small to use the wave functions (JsJ) . This leads 
to the condition to so 3> 3 GeV 2 . At high energy o\ has 
little sensitivity on so which we fix at so = 30 GeV 2 for 
further applications. 

The radiation of each new, n-th gluon can be treated 
as radiation by a color triplet which is an effective quark 
surrounded by n — 1 gluons. It should be resolved by 
the soft interaction with the target to be different from 
the radiation of n — 1 gluons i. e. the radiation cross sec- 
tion is proportional to the difference between the total 
cross sections of the two subsequent Fock states which is 
9C/46o- This can be also proved using a 1/N C expansion 
and the dipole representation of Mueller |l6| . Since the 
radiation of a gluon with a G <§; 1 does not affect the 
impact parameter of the radiating quark, all the quark 
lines in the final state cancel with the same lines in the 
initial state (see the prescription for calculating the ra- 
diative cross section in 0), except for the radiation of 
the n-th gluon. Thus, a n for quark-proton interaction in 
the leading-log approximation reads, 
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Summing up the powers of logarithms in ([!]) we arrive 
at the following expression for the total cross section, 
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and i7q P = (Tq P — 9C/46q. The soft Pomeron intercept, 
ap(0) — 1 + A, and can be evaluated provided that the 
QCD coupling a s is known. 

In Gribov's confinement scenario, chiral symmetry 
breaking occurs when the running coupling a s exceeds 
the critical value a s — a c ~ 0.43 ||. This should hap- 
pen at a distance of the order of the size of a constituent 
quark ~ 0.3 fm. Therefore, this value can be used in (|^). 

One can also calculate the mean (a s ) for nonperturva- 
tive gluon radiation averaging over transverse momenta 
kx of the radiated gluons. The popular way to extend 
the running QCD coupling a s (k 2 ^) down to small kx is a 
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shift of the variable _fc£ =>■ k^ + k 2 , where k 2 « 0.25 GeV 2 
was evaluated in J10| using the dispersive approach to 
calculating higher twist effects in hard reactions Jl8| . 
The nonperturbative interaction of the radiated gluons 
drastically suppresses small transverse momenta, push- 
ing (kj*) to higher values which lowers a s . We use the 
transverse momentum gluon distribution calculated in || 
in the light-cone approach in terms of the universal color 
dipole cross section O], We calculated (a s ) with a sim- 
ple parameterization <j(p) oc 1 — exp(p 2 / . For a rea- 
sonable variation of po = 0.3 — 1 fm the mean coupling 
is in the range (a s ) = 0.38 — 0.43 which is very close the 
the critical value mentioned above 0. Taking the mid 
value (a s ) — 0.4 we get from M), 



A = 0.17 ±0.01 
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This value is about twice as large as the one suggested 
by the data for the energy dependence of total hadronic 
cross sections jL2| . However, the radiative part is a rather 
small fraction of the total cross section (at medium high 
energies). A structure similar to ([?]) with a large A was 
suggested in j2j|] (with quite a different motivation) and 
proved to agree well with the data. 

The factor C in the second term in (0) can also be 
evaluated. We calculated the dipole cross section with 
the gluon effective mass 0.15 GeV (to incorporate con- 
finement) and a s = 0.4 and found C = 2.3 at p = l/£>o- 
Thus, the energy dependent term in (0) is fully deter- 
mined. 

The cross section (Q) apparently violates the Froissart 
bound and one should perform unitarity corrections. In- 
deed, the partial elastic amplitude shows a precocious 
onset of unitarity restrictions at small impact parame- 
ters important even at medium high energies JTs] ] . 

Following jl2|,^0| we assume that the i-dependence of 
the pp elastic amplitude is given by the Dirac electro- 
magnetic formfactor squared. Correspondingly, the mean 
square radius (f 2 h ) evaluated in should be smaller 
than (r 2 ch ). 

For the dipole parameterization of the formfactor the 
partial elastic amplitude which is related via unitarity to 
afP, given by (||), (||), takes the form, 
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where K^(y) is the third order modified Bessel function 
and y = b\j8/B n . The slope parameter grows linearly 
with n due to the random walk of radiated gluons with a 
step 1/6q in the impact parameter plane, B n — 2(r 2 h }/3+ 
n/2bl 

We unitarize the partial amplitude Im jp(s,b) = 
^2 Im7 n (s, b) using the quasi-eikonal model |2l|], 
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where D(s) — 1 = <r s d{s) / a e i(s) is the ratio of the single 
diffractive to elastic cross sections. It is approximately 
equal to 0.25 at the lowest ISR energy and slightly de- 
creases with energy oc s -0 04 |2^j2^]. Note that good 
results can be also achieved with a different unitariza- 
tion scheme similar to one suggested in Q . The details 
will be presented elsewhere. 

In order to calculate the total cross section, cr tot = 
2 J d 2 6Imr(6, s), one needs to fix the energy independent 
term with n = in (|l0|). This can be done comparing 
with the data for a to t at any energy sufficiently high to 
neglect Reggeon contributions. We used the most precise 
data (H at y/s = 546 GeV and fixed a = 39.7 mb. 



The predicted energy dependence of a^ t is shown by 
the dashed curve in Fig. [l] which is in good agreement 
with the data at high energy p(|, but apparently needs 
Reggeon corrections towards low energies. We added a 
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FIG. 1. Data for total pp (full circles) and pp (open circles) 
cross section |^] and the prediction of Eq. ([u)|) for the energy 
dependence of the Pomeron part (dashed curve). The solid 
curves include Reggeon contributions fitted to the data. 



Reggeon term Imffl(s,fe)[l — ImFp(s,6)] screened by 
unitarity corrections, which was fitted independently for 
pp and pp, <7^ p = 17.8 mb/ \J s/sq, = 32.8 mb/ ys/so- 



The fitted Reggeon slope is Br 
where a' R — 0.9 GeV 



r2 r 



2 and R% = 3 GeV~ 



2a' fl ln(s/s ), 



The results are shown by the solid curves of Fig. |l| (pp 
bottom curve and pp upper curve) . 

As soon as the partial amplitude ( |TT| ) is known, we 
are in position to predict the slope of elastic scattering 
at t = 0, B e i(s) = (b 2 )/2, where averaging is weighted 
by the partial amplitude ([ll]). The results exhibit good 
agreement when compared with the pp and pp data ]26| ] 
in Fig. ||. Although the value of the slope essentially de- 
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FIG. 2. Data for the elastic slope |26| and our predictions. 
The upper and bottom curves and open and full circles cor- 
respond to pp and pp respectively. 



pends on our choice of {f 2 h ) in (|l0|), the predicted energy 
dependence, i.e. the effective value a' P is fully defined 
by the parameter bo fixed in ||. Indeed, each radiated 
gluon makes a "step" ~ 1/6q = (0.3 /m) 2 in the impact 
parameter plane leading to the rising energy dependence 
of the clastic slope. Eventually, at very high energies the 
approximation of small gluon clouds breaks down but 
the mean number of gluons in a quark (n) = Aln(s/so) 
remains quite small even in the energy range of collid- 
ers. It is only (n) = 0.7 — 1 at the ISR and reaches 
about two gluons at the Tevatron. Correspondingly, the 
mean square of the quark radius grows from 0.06 fm 2 
to 0.18 fm 2 which is still rather small compared to the 
mean square of the charge radius of the proton. 

Summarizing, the strong nonperturbative interaction 
of radiated gluons substantially shrinks the gluon clouds 
around of valence quarks. This spots are small (~ 
0.3 fm) compared to the hadronic radius, but the gluon 
radiation grows with energy as s A where A = 0.17±0.01. 
Such a steep rise does not contradict the data since this 
fraction of the total cross section is rather small (it con- 
tains a factor X/oq w 1 mb). A large energy independent 
fraction comes from the Born term which corresponds to 
scattering of the valence quark skeleton without gluon 
radiation. A very soft interaction which cannot resolve 
and excite the small spots contributes to this term. It 
cannot be reliably predicted and is fixed by data , while 
the energy dependent term is fully calculated. The re- 
sults are in good agreement with the data for total pp 
and pp cross sections and elastic slopes. 

Note that although we have some room for fine-tuning 



in the the parameters (C, sq, (r 2 h )), the results are 
rather insensitive and the agreement with data is always 
pretty good. We have also tried a different unitarization 
scheme suggested in |24| ] arriving to similar results. 

The details of calculations and further comparison with 
elastic scattering data will be published elsewhere. 
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